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The field induced quantum critical properties of the three dimensional spin- 1/2 anisotropic antifer- 
romagnetic Heisenberg model has been studied. We have investigated the quantum phase transition 
between the spiral order and field induced ferromagnetic order by means of Bose-Einstein conden- 
sation of magnons in terms of a bosonic representation. The effect of in-plane anisotropy on the 
critical properties has been studied via the bosonic model by Green's function approach. We have 
^ , found an analytic expression for the gap exponent in addition to numerical results for the critical 

^ O . magnetic field in terms of anisotropy parameter. The in-plane anisotropy breaks the U(l) symmetry 

' explicitly which changes the universal behavior by a drastic change on the gap exponent. Moreover, 

the critical magnetic field depends strongly on the in-plane anisotropies. The divergence of the 
transverse structure factor at the antiferromagnetic wave vector confirms the onset of the magnetic 
order which scales with the negative value of gap exponent as the magnetic field approaches the 
critical one. The transverse staggered magnetization as an order parameter vanishes with exponent 
/3 = 0.5 when the magnetic field reaches its critical value in low field region. 
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PACS numbers: 75.10.Jm, 75.40.Cx, 75.30.Kz 



I. INTRODUCTION 



Quantum phase transitioni'-S. is an interesting topic for both theoretical and experimental condensed matter research 
activities. This phase transition is found at the zero temperature based on variation of non thermal control parameter 
Q ' such as magnetic field or hole doping. Quantum phase transition (QPT) occurs at the quantum critical point where 
quantum fluctuations destroys a long range order of the model at absolute zero temperature. One of the novel 
species of quantum phase transition is field-induced-magnetic phase transition that can be occurred in insulating 
antiferromagnetic systems^ like transition metal oxides and local spin systems. This kind of QPT has been observed 

■ in the copper halide CS2CUCI4 which is an insulator and each Cu^"*" carries a spin of 1/2. 

t~~«- CS2CUCI4 can be described as a quasi- two-dimensional spin 1/2 antiferromagnet on a triangular lattices {he plane) 

10 weakly coupled along the crystallographic a direction.'* The crystal field effects quench the orbital angular momentum; 
[ however, the anisotropic effect is still significant on the phase transition. The layered crystal structure confines the 
' main superexchange routes to neighboring spins in the he plane. According to the above facts the magnetic properties of 
this material can be described by the antiferromagnetic Heisenberg model. The magnetic field is applied perpendicular 
to the he plane which adds a Zeeman term to the model. For magnetic fields {B) close to the critical field [Be), Zeeman 
term competes with the spin exchange interaction and system enters a field induced ferromagnetic state^i^. The field 
induced ferromagnetic phase {B > Be) has gapped quasi-particles, gapped magnons. The field induced gap vanishes 
at Be when the magnetic field is reduced and the magnetic ordering for transverse component of spins sets up. This 
latter state (for B < Be) is named spiral long range order. In this work we have studied the mentioned quantum 

■ phase transition based on Bose-Einstein condensation of magnons via a bosonic gas model^. Bloch- applied the 
" " ' Bose-Einstein quantum statistics to the excitations in solids which gives the basic notion to relate Bose-Einstein of 

magnons to the magnetic ordering in the original spin model^i^. 

The isotropic Heisenberg model with longitudinal applied magnetic field has been studied by theoretical and numer- 
ical methods. The isotropic model on cubic lattice has been investigated by numerical quantum Monte-Carlo method 
at finite temperature which gives the phase boundary between spiral order and induced ferromagnetic state*°. A 
theoretical approach based on Bose condensation of magnons for the isotropic model on triangular lattice has been 
studied in Ref. [§1 . Furtheremore, the experimental data for specific heat indicate that the A like anomaly peak appears 
in the behavior of specific heat versus temperature for B < Bc{T ~ 0)^. 

In the general case, the spin model Hamiltonian can include spatial anisotropies in the exchange coupling between 
nearest neighbor spins. This property is related to the existence of easy axes magnetization due to crystalline 
electric field and spin-orbit coupling. Dzyaloshinskii-Moriya (DM) interaction with a DM vector in a specific direction 
establishes easy-plane spin anisotropy in the Cs2CUCl4*^. Anisotropy due to DM interaction violates SU(2) symmetry 
of the isotropic Hamiltonian although U(l) symmetry corresponding to spin rotation around DM vector is still present. 
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However, the spin-orbit coupling may induce the anisotropy in the be easy-plane which reduces the U(l) to Z(2) 
symmetry. Therefore, Goldstone theorem^^ can not be applied in such cases, because it is applied for Hamiltonian 
with a continuous symmetry. In other words, the excitations of this model are not Goldstone modes. 

In this paper, we intend to find the effect of in plane anisotropy on the critical point and the transverse spin 
structure factor close to the field induced QPT. We have considered the fully anisotropic spin 1/2 Heisenberg model 
in the presence of a longitudinal field on a cubic lattice. We anticipate that the general behavior for cubic lattice is also 
valid for the triangular one. Moreover, the study on cubic lattice reduces the complexity of calculations which will be 
the route to investigate the triangular case. In addition the results on cubic lattice can be applied to the field-induced 
magnetic phase transition of TlCuCla . We have implemented the hard core boson transformation for spin operators 
which gives the excitation spectrum in terms of many body calculations for bosonic gasi^. We have used Brueckner 
approach^i to find the bosonic self energy to get the magnon dispersion relation. The quantum critical point is 
approached where the magnon spectrum becomes gapless. We have found an analytic expression for the gap exponent 
in terms of the anisotropy parameter. Our results show that a small amount of in-plane anisotropy changes the gap 
exponent drastically which is the witness for the change in universal behavior. We have also found the dependence of 
critical magnetic field on the anisotropy parameter which is also justified by the divergence of the transverse structure 
factor at the antiferromagnetic wave vector. The divergence of the in-plane magnetic susceptibility obeys an algebraic 
power law with an exponent equals to the negative of gap exponent as the magnetic field approaches the critical 
one. Moreover, the vanishing of the staggered magnetization is given be the exponent /? = 0.5 in the mean field 
approximation. 



The most general effective Hamiltonian to describe magnetic insulating matter due to exchange interaction between 
the spin of localized electrons can be written by 



where g ~ 2.2, fis is the Bohr magneton and B is the magnetic field. The localized spins are located on the cubic 
lattice structure with nearest neighbor exchange interaction. The effect of spin orbit coupling is generally entered to 
the fully anisotropic exchange couplings J,". The exchange anisotropy is defined by parameter v with the following 
relations 



where the scale of energy J is set to one. The above type of exchange interaction implies an anisotropy in both in 
plane and axial directions where the symmetry has be reduced to Z(2). 

III. HARD CORE REPRESENTATION OF SPIN HAMILTONIAN AND BOSONIC GREEN'S 

FUNCTIONS 

As mentioned in the introduction we intend to describe the field induced QPT in terms of Bose-Einstein condensation 
of magnons. Our approach is similar to what we have implemented in RefllSl to study the quantum critical properties 
of the Kondo-necklace model. In the first step the spin Hamiltonian is mapped to a bosonic model. This is done 
by the hard core boson representation given by: — > ai, — > a| and S'f = 1/2 — ajai, where Oj and a] are 
boson annihilation and creation operators, respectively. The SU(2) algebra of spin operators is retrived from bosonic 
algebra of and operators with the hard core constraint, i.e only one boson can occupy a single site of lattice. 
The constraint is added to the Hamiltonian by an on-site infinite repulsion among bosonic particles. The resulting 
Hamiltonian in terms of bilinear {Hbu) and interacting {Hint) parts is given by 



II. ANISOTROPIC SPIN HAMILTONIAN 




(1) 



(ij) a=x,y,z 



r = j(i + u) , jy = j(i -v) , = J, 



(2) 




(3) 




(4) 
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The bilinear Hamiltonian in the Fourier space representation is 



(5) 



^k 



(6) 



It is defined Jk — J'^a=x y % cos(fcc() and Jq = 3 J. The wave vectors fca are considered in the first Brillouin zone. 
The effect of hard core repulsion part — ^ oo) of the interacting Hamiltonian in Eq.(j4|) is dominant compared with 
the second quartic term. Thus, it is sufficient to take into account the effect of hard core repulsion on the magnon 
spectrum and neglect the second quartic term. The interacting part of Hamiltonian in terms of Fourier transformation 
of bosonic operators is given by 



k,k' ,q 



(7) 



The bilinear Hamiltonian is simply diagonalized by the unitary Bogoliuobov transformation to the new bosonic quasi 
particle operators Ofk and aj^, 15 1, which is given by 



y-bii = y^^k (alak + 1/2) , 



uji = Ai-B? 



k7 



and the Bogoliubov coefficients are 



- - + — • 



(8) 



(9) 



Although the bilinear part of Hamiltonian is diagonal in the new bosonic (ak, cejj.) representation, to avoid the 
complexity of the calculations for the hard core repulsion term the Green's functions will be calculated in the original 
boson operators (afc,a|,). In the original boson representation, Hbu includes the pairing term between magnons 
which requires both anomalous and normal Green's functions to be considered. More explanations of the detailed 
calculations can be found in Ref, 15. Finally, the self-energy is expanded in the low energy limit which gives the single 
particle part of Green's function (Gf^), 



ZkUl 



ZkVi 



(10) 



where the renormalized triplet spectrum (fife), the renormalized single particle weight constants {Zk) and renormalized 
Bogoliuobov coefficients ([/fe, T4) are given by 

= Zk^[Ak + S„(fc,0))]2 - [Bk + Sa(fc,0)]2), 



'-•It' 



^ ^1 ^ Zfc[^fc + S„(fc,0)] 
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(11) 



The renormalized weight constant is the residue of the single particle pole in the Green's function. In the next step 
we will take into account the effect of hard core repulsion on the magnon spectrum. 

IV. EFFECT OF HARD CORE REPULSION ON THE MAGNON SPECTRUM 



The density of the magnons is obtained from the normal Green's functions 

rii = {a\ai) ^ -j^^vl, 

k 



(12) 



4 



where N is the number of the spins in the cubic lattice. In the vicinity of the critical field (-Bj?) and at the zero 
temperature the density of excited magnons is negligible^-. Since the Hamiltonian Hint in Eq.([7]) is short ranged 
and U is large, the Brueckner approach (ladder diagram summation) -'^^i^^ can be applied for the low density limit of 
magnons. The interacting normal Green's function is obtained by imposing the hard core boson repulsion, U — >■ oo. 
Firstly, the scattering amplitude (t-matrix) r(fci, fc2; ^3, ^4) of magnons is introduced where ki = (k, (fco))i. The 
basic approximation made in the derivation of T(K) is that we neglect all anomalous scattering vertices, which are 
presented in the theory due to the existence of anomalous Green's functions. According to the Feynman rules^ in 
momentum space at zero temperature, the scattering amplitude is calculated (see Fig.l of Ref,15). By replacing the 
noninteracting normal Green's function in the Bethe-Salpeter equation and taking the limit U — > 00 we obtain the 
scattering matrix in the form 

r(K,.) = -(^ / d^Q "^""-^ ^^^y\ (13) 

V(27r)'^J W-CJQ-WK-Q LJ + UJQ + UJ-K-Q / 

According to Fig. 2 of RefflBI and after some calculations the normal self-energy is obtained in the following form 

S^(k,^) = |^i;^r(p + k,c.-c.p). (14) 
p 

In the dilute gas approximation there are other diagrams which are formally at most linear in the density of magnons. 
However, the magnon densities are very small and the contributions of such terms are numerically smaller than 
Eq. (IT4|) . We should also consider the anomalous self-energy related to Hu which exists in the vertex function. The 
anomalous self-energy has a vanishing contribution. 



V. THE GAP EXPONENT 



Close to the quantum critical point, the excitation gap (A) in the field-induced ferromagnetic phase vanishes 
according to the following power law behavior 

A^|B-Be|^ (15) 

where Be is the critical magnetic field and (f> is the gap exponent which is related to universality class of the quantum 
critical point. The quantum critical point corresponds to the vanishing of magnon spectrum at the antiferromagnetic 
wave vector Qaf = {t^, tt, tt). The magnon spectrum close to the antiferromagnetic wave vector (Qaf) is approximated 

by 



where c is the spin wave velocity. The spin wave velocity (c) is obtained numerically from the excitation spectrum. 
In the first step, we calculate the variation of the self-energy related to Hu which is given by 

S^'^iQAp) = ^^'5i^^r(A; + Q^f.,-c^,) + ^^t|^r(A; + Q^^.,-u;fe). (17) 

k k 

The main contribution to the first integral in Eq. (jl7p comes from the small momenta q ~ A/c ^ 1 where q = k — QAF 
since 

= + ^^^^[-]) - - ,Jl'^L. - (18) 

2\ UJq Wfc / 2(A^ + c'^g^)'^/^ 

Taking into account the first correction to the magnon density, the vertex function can be written for small q (see 
Ref.tlTl]) 

r{q,-u;,)^n[l+^^^^\nq], (19) 

where Fq = r'^{k = 0) and the value of quantity X at the critical point is defined by X'^. The substitution of Eq. p^ 
in Eq. pT]) and replacing g ~ A/J in the first integral of Ea. (IT7|) . we find that 

SJ:^in,n) = -%^rS[l + ^^^ln^]+T'n,SB, (20) 
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where V = ^'"^'^'g"''"' and nh{= the density of magnons at the critical point. Let us define the following 



SB "^^^ '-o\— N /L^i\^i 

expressions 



A = 



AF- 



(J = T'nb. (21) 
After some calculations we finally get the following relation 

= (^Al£±^ (i _ ^^iJ-^) . (22) 
A V 47rc2 J J 

The gap exponent is obtained upon replacing A = jiJi?!*^ in the above equation which finally is equal to 

* = 5--|^. (23) 

We have presented the numerical results of gap exponent in terms of anisotropic parameters in section- TVIII 

VI. STAGGERED MAGNETIZATION FOR B < Be 

At the quantum critical point the magnons condensate at q = Qaf which is the onset of the long range antifer- 
romagnetic (AF) order in the model. The system is represented by the AF ordered state as far as B < Be- In the 
AF phase we have applied Hartree-Fock-Popov mean field approach^i^^ by taking into account the condensation of 
magnons in the interacting Hamiltonian at wave vector Qaf- The effective interparticle interaction is defined by the 
following Hamiltonian 

Heff = V%Qaf) ' (24) 

where, T'^{Qaf) is the interaction parameter at the critical point [B = Be). Below the critical point {B < Be), 
the AF order parameter becomes nonzero and it can be interpreted by the nonzero mean field value of the creation 
operator of magnons at Qaf- Let us define {aq^^p) — Nc = (ag^^) for the condensate phase where Nc is the number 
of condensed magnons. The staggered magnetization in the x-y plane which represents the long range AF order is 
denoted by m± = nix + inT-y — gfJ-ByNn^ where N is the total number of sites and ric is the condensed magnon 
density. The effective Hamiltonian in Eq.(|24|) can be written in the following form where the contribution from the 
condensate phase has been denoted by Hjj, 

Hu = + + + H^, 

r%QAF)N^ 



2N 



2 ^''{Qaf)Nc sr^ \l , t t \ 1 o t 

q 

Hu = 2^ [alak+qa^q + h.c. j , 

k.q 

rj ^"(Qaf) t t ^OK\ 



k,q,k' 



In the above equations, ^ implies that the terms with creation and annihilation operators at the antiferromagnetic 
wave vector {Qaf) are excluded. In a mean field approximation the contribution form is zero since it contains 
linear terms of boson operators. Taking into account the hard core repulsion which avoids the pairing of magnons 
and considering all other contractions the mean field representation of is 

H^^2{l~n,)T'{QAF)Y,alak. (26) 
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After adding the non- interacting part, Eq.([5]), to the mean field (MF) interacting one the Hamiltonian is given by the 
following equation plus a constant term which has been omitted here, 

H'^^ = j^iAu + 2T^iQ^,))aia, + ^ ^^I^tll^^^^a-, + 4aL,). (27) 

k k 

The mean field Hamiltonian is diagonalized by the unitary Bogoliubov transformation. 



= ViAk + 2r-{QAF))^ - {vJk + r=(gAF)nc)2, 

flfc = dk(f)k - fk<P^_k' (28) 

where fl^ gives the excitation spectrum of the new bosonic quasi partiles defined by the creation operator (/)^ and d^, 
fk are Bogoliubov coefficients 



'^^^V %h +2' ■^^^ = V 2- (2^) 

The condensation of magnons at the AF wave vector implies that excitation spectrum must be gapless at Qaf which 
gives the following relation 

Jk=QAF - Jo + msB + 2T%af) = -yJk=Q^^ - V'iQAFW (30) 
Therefore, the transverse order parameter for B < Be is given by 



^ / - Jk=QAF + JS ~ g^B-^ - 2r^(Q^f ) - i^Jq^^ 
m± = guBV^^ gfJ-Bd V^iQAF) 

The critical field is touched where ric = in Eg. ([50)) which leads to the following expression, 

^ ^ Jq - ^Qaf - 2r^(QAF) - i^JQAF ,32) 

Therefore the transverse staggered magnetization has the following expression in the mean field approximation 



^ gfisiBc-B) , ^ 

The scaling behavior of transverse order parameter close to the critical field is characterized by the exponent /3 via 
m_L ~ \Bc — B\^ which is /3 = 0.5, in the mean field approximation. 



VII. SUMMARY, RESULTS AND DISCUSSIONS 



In this article we have studied the effect of in-plane anisotropy on the quantum critical properties of the spin 1/2 
Heisenberg model in the presence of a longitudinal field {B) on a cubic lattice. The in-plane anisotropy breaks the 
U(l) symmetry of the model (around the direction of the magnetic field) and changes the quantum critical properties 
of the model which is discussed in this section. Moreover, we have analyzed the field-induced phase transition in this 
model in terms of Bose-Einstein condensation of magnons. 

The original spin model has been represented by a bosonic model in the presence of hard core repulsion to avoid 
double occupation of bosons at each lattice site which preserves the SU(2) algebra of the spin model. In the limit of 
B/J — > 00, the ground state is a field induced ferromagnetic state and a finite energy gap exists to the lowest excited 
state which is called the magnon spectrum. The decrease of magnetic field lowers the excitation gap which eventually 
vanishes at the critical magnetic field (Be). This point corresponds to the condensation of magnons which is the 
onset of long range antiferromagnetic order of the spin model. We have implemented the Green's function approach 
to obtain the effect of interaction on the diagonal part of the bosonic Hamiltonian using Brueckner formalism close 
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FIG. 1: (Color online) The energy gap versus the magnetic field for various anisotropy parameters. The change of the critical 
magnetic field (where the gap vanishes) for various anisotropics is remarkable. 

to the quantum critical point where the magnon density is small. The magnon spectrum have been calculated from 
Eas. (fT0l[Tn[T4|) selfconsistently. The procedure is started with an initial guess for Z^, T,n{k, 0) and T,a{k, 0), then using 
Eg. dlip we find the renornialized excitation energy and the renormalized Bogoliuobov coefficients. The procedure is 
repeated until convergence is reached. Using the final values for energy gap, renormalization constants and Bogoliubov 
coefficients, we have obtained the quantum critical point for different anisotropy parameters in TableU Our data 
shows that a small amount of anisotropy has a considerable change on the critical magnetic field. We have also plotted 
the magnon gap versus the magnetic field in Fig |1] for different values of anisotropy. It is obvious from Fig |1] that the 
magnon gap vanishes as the magnetic field approaches the critical value for B ^ B^- 

Moreover, the scaling behavior of gap close to Be which is characterized by the gap exponent (0) defined in Eq.([T5]) 
depends on the anisotropy parameter. We have presented the gap exponent (0) for different anisotropics in TablelH 
The dependence of on z/ shows that the in-plane anisotropy changes the universal behavior of the model. A drastic 
change of from 0.40 for i/ = to 0.2 for v = 0.1 manifests the change of universality class due to explicit breaking 
of symmetry by anisotropy. hi v = the model has U(l) symmetry while for v ^ Q the symmetry breaks to Z2. 
Although the calculated gap exponent for different nonzero anisotropics may change slightly, the amount of changes 
are in the order of error bar which mainly comes from the error bar inherited in the value of Bc{i'), nonzero values of 
V do not present different universal behavior. 



TABLE I: The critical magnetic field (-Be) and gap exponent (0) for different values of anisotropics. The error bar for all data 
is ±0.05. 



V 


0.0 


0.1 


0.2 


0.3 


0.4 


Be (critical field) 


12.26 


13.03 


13.88 


14.72 


15.53 


(gap exponent) 


0.40 


0.2 


0.2 


0.2 


0.2 



The long range ordering can be deduced from the behavior of spin susceptibility. The magnetic ordering occurs for 
transverse spin components; therefore, static spin structure factor for tranverse component diverges at the antiferro- 
magnetic wave vector in the quantum critical point. The x-component static spin structure factor at momentum q is 
defined by 




X"^(g) = (s^(g)s-(-g)), 



(34) 
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FIG. 2: (Color online) The (x-component) transverse structure factor at the antiferromagnetic wave vector versus the magnetic 
field for difi'erent anisotropics. The divergence at the critical magnetic field justifies the onset of magnetic order. 



which is given by the following expression, 



Close to the quantum critical point where the magnon spectrum vanishes the dominant contribution is obtained to 
be 



X^^(QAF)«7r^^,Ji3-Bc|-^ (36) 

which shows that the divergens of magnetic susceptibility follows a scaling relation with exponent, (p. The numerical 
results for x^^{Qaf) versus the magnetic field have been plotted in Figl2] This plot confirms the presence of antifer- 
romagnetic order at the critical magnetic field. The divergence of the static structure factor happens at the different 
critical fields for various anisotropics {v) which justifies the previous results on energy gap. 

Both results on the energy gap and transverse structure factor confirms that a small amount of in-plane anisotropy 
changes the critical magnetic field considerably. Moreover, the explicit breaking of symmetry due to in-plane 
anisotropy shows up in the gap exponent. Our results on the scaling behavior of magnetic order parameter (/?) 
is limited to the mean field approximation which does not show its dependence on the anisotropy. However, we expect 
that the dependence of (3 on anisotropy {v) should appear if the calculation goes beyond the mean field approach. 
As far as the model has two control parameter v and B, the universal behavior should be fixed by two exponents <j) 
and j3. In other words, any other exponent for a scaling behavior close to critical field at zero temperature can be 
expressed in terms of the obtained exponents ((/>, (3). 
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